UNCLASSIFIED 


ADA409975 


The  Ideal  Constant  Volume  Limit  of  Pulsed  Propulsion 


AIR  FORCE  RESEARCH  LAB  EDWARDS  AFB  CA 

25  OCT  1999 

Distribution  Statement  A: 

Approved  for  public  release.  Distribution  is  unlimited. 


UNCLASSIFIED 


UNCLASSIFIED 


Redistribution  Of  DTIC-Supplied  Information  Notice 

As  a  condition  for  obtaining  DTIC  services,  all  information  received  from  DTIC  that 
is  not  clearly  marked  for  public  release  will  be  used  only  to  bid  or  perform  work 
under  a  U.S.  Government  contract  or  grant  or  for  purposes  specifically  authorized 
by  the  U.S.  Government  agency  that  sponsored  the  access.  Furthermore,  the 
information  will  not  be  published  for  profit  or  in  any  manner  offered  for  sale. 

Reproduction  Quality  Notice 

We  use  state-of-the-art,  high-speed  document  scanning  and  reproduction  equipment. 
In  addition,  we  employ  stringent  quality  control  techniques  at  each  stage  of  the 
scanning  and  reproduction  process  to  ensure  that  our  document  reproduction  is  as 
true  to  the  original  as  current  scanning  and  reproduction  technology  allows. 
However,  the  following  original  document  conditions  may  adversely  affect 
Computer  Output  Microfiche  (COM)  and/or  print  reproduction: 

•  Pages  smaller  or  larger  than  8.5  inches  x  11.0  inches. 

•  Pages  with  background  color  or  light  colored  printing. 

•  Pages  with  smaller  than  8  point  type  or  poor  printing. 

•  Pages  with  continuous  tone  material  or  color  photographs. 

•  Very  old  material  printed  on  poor  quality  or  deteriorating  paper. 


If  you  are  dissatisfied  with  the  reproduction  quality  of  any  document  that  we 
provide,  particularly  those  not  exhibiting  any  of  the  above  conditions,  please  feel 
free  to  contact  our  Directorate  of  User  Services  at  (703)  767-9066/9068  or 
DSN  427-9066/9068  for  refund  or  replacement. 


Do  Not  Return  This  Document  To  DTIC 


UNCLASSIFIED 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  tor  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  this  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis 
Highway,  Suite  ,1204,  Arlington,  VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a 
collection  of  information  if  it  does  not  display  a  currently  valid  OMB  control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 

1.  REPORT  DATE  (DfrMM-YYYY)  27REPORT  TYPE  ~~  3^ DATES  COVERED  fF/mi  ^To) 

_  Technical  Papers _ 

4.  TITLE  AND  SUBTITLE  5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 


5c.  PROGRAM  ELEMENT  NUMBER 


6.  AUTHOR(S) 


5d.  PROJECT  NUMBER 

.3  05 1 

5e.  TASK  NUMBER 

gEl IB 

Sf.  WORK  UNIT  NUMBER 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Air  Force  Research  Laboratory  (AFMC) 

AFRL/PRS 
5  Pollux  Drive 

Edwards  AFB  CA  93524-7048 

9.  SPONSORING /MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


Air  Force  Research  Laboratory  (AFMC) 
AFRL/PRS 
5  Pollux  Drive 

Edwards  AFB  CA  93524-7048 

12.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 


8.  PERFORMING  ORGANIZATION 
REPORT 


10.  SPONSOR/MONITOR’S 
ACRONYM(S) 


11.  SPONSOR/MONITOR’S 
NUMBER(S) 


Approved  for  public  release;  distribution  unlimited. 


13.  SUPPLEMENTARY  NOTES 


14.  ABSTRACT 


20030123  036 


15.  SUBJECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION 

18.  NUMBER 

OF ABSTRACT 

OF  PAGES 

a.  REPORT 

b.  ABSTRACT 

c.  THIS  PAGE 

(O) 

Unclassified 

Unclassified 

Unclassified 

...  ■  UV- . 

>p  ay 


PERSON 

Leilani  Richardson 
19b.  TELEPHONE  NUMBER 

(include  area  code) 

(661)275-5015 

Standard  Form  298  (Rev.  8-98) 
Prescribed  by  ANSI  Std.  239.19 


MEMORANDUM  FOR  PR  (In-House  Publication) 

FROM:  PROI  (TI)  (STINFO)  25  October  1999 

SUBJECT:  Authorization  for  Release  of  Technical  Information,  Control  Number:  AFRL-PR-ED-TP-1999-0200 
Talley,  D.,  “The  Ideal  Constant  Volume  Limit  of  Pulsed  Propulsion” 


36th  AIAA/ASME/SAE/ASEE  Joint  Propulsion  Conference  and  Exhibit  (Statement  A) 


SOME  ANALYTICAL  SOLUTIONS  FOR  THE 
IDEAL  CONSTANT  VOLUME  LIMIT  OF  PULSED  PROPULSION 


Doug  Talley 

Propulsion  Directorate,  Air  Force  Research  Laboratory 
Edwards  AFB,  CA 

ABSTRACT 

A  constant  volume  limit  of  pulsed  propulsion  is  formulated  which  applies  when  blow  down 
times  are  much  longer  than  characteristic  wave  transit  times  in  the  combustion  chamber.  Under 
this  limit,  the  combustion  chamber  is  approximated  as  being  time  varying  but  spatially  uniform, 
and  the  nozzle  is  approximated  as  being  one  dimensional  but  quasi-steady.  Some  analytical 
solutions  for  this  limit  with  fixed  expansion  ratio  nozzles  are  explored  for  the  isentropic  blow 
down  of  a  constant  y  ideal  gas.  The  results  are  compared  with  a  variable  expansion  ratio  case 
where  the  expansion  ratio  is  continuously  varied  to  match  the  pressure  ratio  during  blow  down. 
The  major  conclusions  are  that  constant  volume  devices  should  optimize  at  approximately  the 
same  mixture  ratios  as  constant  pressure  devices,  and  that  using  a.  fixed  expansion  ratio  results  in 
only  a  modest  impulse  penalty,  not  exceeding  3%  for  the  cases  examined,  compared  to  using  a 
variable  expansion  ratio,  as  long  as  the  fixed  expansion  ratio  has  been  optimized  to  produce  the 
maximum  possible  impulse  for  die  blow  down. 

1.  INTRODUCTION 

Recent  interest  in  pulsed  detonation  propulsion  has  spawned  a  number  of  attempts  to  model  the 
system  performance,  the  results  from  which  have  so  far  varied  widely  [1].  Constant  volume  (CV) 
combustion,  as  refered  to  here,  is  a  limiting  case  for  pulsed  combustion  cycles  which  is 
approached  when  blow  down  times  become  much  longer  than  characteristic  wave  transit  times  in 
the  combustion  chamber  (large  combustion  chambers  and  small  nozzles).  In  this  limit,  the 
combustion  chamber  can  be  approximated  as  being  time  varying  but  spatially  uniform,  and  the 
nozzle  flow  can  be  approximated  as  being  one  dimensional  but  quasi-steady.  Compared  with 
how  a  real  device  is  likely  to  operate,  the  CV  limit  probably  underpredicts  the  thrust 
(unrealistically  small  nozzles),  but  the  specific  impulse  is  probably  close  or  equal  to  an  upper 
bound.  Thus  the  CV  limit  is  a  useful  reference  case.  Some  analytical  solutions  for  the  CV  limit 
are  explored  here  for  the  isentropic  blow  down  of  a  constant  y  ideal  gas  at  fixed  expansion  ratio. 
The  results  are  compared  with  the  case  where  the  expansion  ratio  is  continuously  varied  to  match 
the  pressure  ratio  at  all  times  during  the  blow  down. 

2.  FORMULATION 

2.1  General 

The  impulse  produced  by  the  unsteady  blow  down  of  the  combustion  gases  in  the  CV  limit  is  the 
integral  of  the  thrust  F  =  mve  +  {Pe  -  P„)Ae  over  time,  where  m  is  the  mass  flow  rate,  Pm  is  the 
ambient  pressure,  and  ve  ,  Pe,  and  Ae  are  the  velocity,  pressure,  and  area,  respectively,  at  the 
exit  plane.  Under  the  transformation  dt  -  {dt  /  dp)dp  -  -(V  /  m)dp  and  the  further 


transformation  r  =  p/pQ,  where  V  is  the  combustion  chamber  volume  and  p0  is  the  initial 

chamber  density  after  combustion  but  before  blow  down  at  timer  =  0,  the  total  impulse  I  and 
blow  down  time  t  may  be  expressed  as 


m  =  J/(r)A  =  p0Vj(ve  +  (Pt  -  PJA,  hh)dr , 
t  =  ^dt  =  p0V^  dr/ m 


2.2  Fixed  Nozzles 


For  isentropic  flow,  the  instantaneous  exit  velocity  is  v,  =  j2cp(T  -Tt) ,  where  T  and  Te  are 

the  instantaneous  chamber  and  exit  temperatures,  respectively.  Pressures  and  densities  in  both 
the  chamber  and  the  nozzle  will  be  related  in  a  simple  fashion  to  temperatures  according  to 
T IT0  =  (p /  p0)r"‘  =  (PI P0)(r-!)/r ,  where  the  subscript  “0"  denotes  conditions  in  the  chamber  at 

time  t  =  0 .  With  c0  =  ^yRT0  and  cp=yR  /(y  - 1) ,  the  exit  velocity  can  be  expressed  as 


v,=V' 


where  y  is  the  ratio  of  specific  heats  and  re  =  pj  p  is  the  ratio  of  the  density  at  the  exit  plane  to 
the  density  in  the  chamber.  For  isentropic  nozzle  flow,  the  exit  density  ratio  is  related  to  the 
expansion  ratio  e  s  Aef  A',  where  A*  is  the  throat  area,  by  ([2],  eq.  5.3) 


,  _  g(r) 


f  y — 1 

s(y)=P~ 


(l/2)(y+])/(r-l> 


Equation  (4)  has  two  solutions  corresponding  to  subsonic  and  supersonic  flow  at  the  exit  plane, 
and  these  solutions  will  depend  only  on  £  and  y .  For  supersonic  flow,  the  density  ratio  rt  will 
therefore  remain  constant  during  blow  down  as  long  as  compression  waves  do  not  enter  the  noz¬ 
zle.  The  exit  velocity  term  in  (1)  will  therefore  be  analytically  integrable  as  long  as  this  is  true. 

For  the  supersonic  solution,  the  choked  throat  conditions  will  be  functions  only  of  r  and  y .  For  a 
throat  density  and  velocity  given  by  ([2],  eq.2.35)  p*  =  rp0[2/(y  +  l)]I/(r_!)  and 

v*  =^jyRT'  =c0(2/(y  +  l))1/2r7'’,  and  for?,  =P0r/rr,  the  instantaneous  mass  flow  and  thrust 
become 


m(r)  =  p*AV  =  p0c0A*g(y ) 


F(r)  =  P0A*g(7)  rry 


7-1 v 


l  -r!A  + 


— -  (~Yrr 


where  (4)  has  also  been  used.  Thus  all  terms  in  (1)  will  integrable  analytically. 

Define  <pG  =  PJ  P„  to  be  the  ratio  of  the  initial  pressure  at  time  t  =  0  to  the  ambient  pressure, 
and  define  i  =  I  /  p0cQV  and  x  =  tc0A *  IV  to  be  the  normalized  impulse  and  the  normalized  time, 
respectively.  Then  integrating  (1)  and  (2)  gives 

o  \U2f  o  's 
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b(r), 


(8) 


where 


.  a(r)  =  1  -  rinl)n  b(r)s(l/ rf~X)'2  -1.  (9) 

The  quantity  i/t  will  be  related  to  the  average  thrust  F  -I/t.  Noting  that  c2  =yRT0  and 
P0  =  p0RT0 ,  then  combining  constants  after  dividing  leads  to 

yi/x-F  JP0A* &cF,  (10) 

where  cF  is  the  average  thrust  coefficient  for  the  blow  down. 

2.3  Variable  Nozzles 

Variable  nozzles  can  be  envisioned  where  the  expansion  ratio  is  continuously  adjusted  to  match 
the  pressure  ratio  at  all  times  to  produce  the  maximum  possible  impulse.  With  pressures  always 

matched,  only  the  exit  velocity  term  v4  =  ^2cp(T  -Te)  in  (1)  will  be  of  concern.  For  fixed  noz¬ 
zles,  it  was  found  that  the  term  TJT  -  rj~x  was  a  constant,  but  here  the  exit  temperature  is  fixed 
by  the  exit  conditions  Pe  =P„ ,  pe  =  po(l/0o)1/r,  and  Te  =Tfl(l/0o)(r'1Vr .  Therefore  TJT  will 
not  be  constant  in  this  case.  If  the  exit'  area  is  assumed  to  be  varied  to  match  pressures  at  all 
times,  as  might  approximately  be  the  case  with  an  automatically  compensating  nozzle  such  as  a 
plug  nozzle,  the  mass  flow  and  blow  down  time  based  on  a  constant  throat  area  are  still  given  by 
(5)  and  (8),  respectively.  The  normalized  impulse  becomes 

/  -  \l/2 


t(r)  = 


y-i 


f1/rH  -(1  TiTJ^dr. 


(11) 


Eq.  (11)  has  no  convenient  analytical  solution,  but  may  easily  be  integrated  numerically. 


3.  RESULTS  AND  DISCUSSION 


3.1  Effect  of  Thermochemistry 

The  dimensionless  impulse  and  blow  down  time  (7)-(9)  depend  explicitly  on  r  and  implicitly 
only  on  y,  r,  (or  £),  and  <p0.  Of  these,  only  y  is  influenced  by  the  thermochemistry,  but  only 

weakly.  The  major  influence  of  thermochemistry  comes  through  the  initial  speed  of  sound  c0 
used  to  normalize  the  impulse.  Thus  the  specific  impulse  Js  =  ic0  is  maximized  when  c0  is 
maximized.  This  in  turn  implies  that  the  optimum  thermochemistry  is  that  which  maximizes  the 
initial  combustion  temperature  T0  and  minimizes  the  molecular  weight.  The  same  general  guid¬ 
ance  is  known  to  also  apply  to  constant  pressure  devices.  Thus  constant  volume  devices  should 
optimize  at  approximately  the  same  mixture  ratios  as  constant  pressure  devices.  Also,  like  con¬ 
stant  pressure  devices,  the  thrust  is  maximized  by  maximizing  the  initial  chamber  pressure  and 
the  throat  area,  as  implied  by  (10). 

3.2  Fixed  Nozzles 

Much  can  be  understood  about  the  blow  down  of  fixed  nozzles  by  examining  conventional 
steady  state  thrust  coefficient  curves  such  as  can  be  found  in  standard  texts  [3].  The  thrust  coeffi¬ 
cient  is  defined  in  (10),  where  P0  in  the  steady  state  case  is  interpreted  to  be  the  steady  chamber 
pressure  of  a  constant  pressure  device.  A  set  of  these  curves  is  given  in  Fig.  1  for  y  =  1.2 .  Di¬ 
mensionless  isobars  (curves  of  constant  <p0)  initially  increase  with  e ,  reach  a  maximum,  then 
decrease  to  a  minimum  where  a  shock  enters  the  nozzle  and  the  above  formulation  is  no  longer 
valid.  The  curve  for  0O  =  »  reaches  a  maximum  only  at  e  =  <» .  A  curve  connecting  the  maxima 
indicates  the  expansion  ratio  producing  the  maximum  thrust  for  a  given  isobar. 

The  blow  down  of  a  fixed  nozzle  CV  device  proceeds  along  a  vertical  path  of  constant  e  be¬ 
tween  two  isobars.  The  thrust  produced  will  be  some  average  between  the  two  isobars.  Picking 
any  two  isobars  in  Fig.  1,  say  between  <p0  =  50  and  <j)0  =  20,  and  following  the  vertical  line  be¬ 
tween  them  for  various  e ,  it  can  be  envisioned  that  the  average  thrust  will  reach  a  maximum 
neither  at  large  e  nor  at  e  =1,  but  at  some  optimum  e  in  between.  However,  the  average  thrust 
cannot  be  calculated  directly  from  Fig.  1,  because  the  thrust  coefficient  is  proportional  to  the 
thrust  divided  by  the  chamber  pressure,  not  the  thrust  alone.  A  blow  down  from  0O  -  50  to 

<j>0  =  20  is  replotted  in  Fig.  2,  where  the  curve  for  0O  =  20  is  normalized  by  the  same  pressure  as 
for  0o=5O,  making  all  curves  proportional  to  the  thrust.  This  is  accomplished  by  multiplying 
the  thrust  coefficient  of  Fig  1.  by  0.4  for  0O  =  20.  The  curve  labeled  “b.d.”  gives  the  average 
thrust  coefficient  for  the  blow  down  and  is  seen  to  be  an  average  of  the  curves  for  0O  =  50  and 
0O  =  20 ,  as  expected. 

Average  thrust  coefficients  are  plotted  as  functions  of  0O  and  e  for  blow  downs  to  r  =  0.75 , 0.5, 
and  0.25  in  Figs.  3, 4,  and  5,  respectively.  Dimensionless  blow  down  times  are  plotted  as  a  func- 


tion  of  r  for  two  values  of  yin  Fig.  6.  As  can  be  seen  from  (8)  and  Fig.  6,  blow  down  times  are 
independent  of  <j)0,  and  depend  only  weakly  on  y.  Because  ris  fixed  and  independent  of  0O  for  a 

fixed  r,  the  curves  in  Figs.  3, 4,  and  5  will  also  be  proportional  to  the  dimensionless  impulse  i,  as 
can  be  shown  from  (10).  However,' the  constant  of  proportionality  will  be  different  for  each  r, 
increasing  as  r  decreases,  because  t  increases  as  r  decreases.  The  constant  of  proportionality  is 
given  in  the  figures.  The  overall  trend  in  these  figures  is  that  the  average  thrust  coefficient  de¬ 
creases  as  r  decreases.  This  is  due  to  T  increasing  more  rapidly  than  i  as  r  decreases,  because  i 
also  increases  as  r  decreases. 

The  expansion  ratio  which  produces  the  maximum  impulse  can  in  principle  be  found  by  taking 
the  first  derivative  of  (7)  with  respect  to  re  and  setting  the  result  equal  to  zero.  The  resulting  ex¬ 
pression  cannot  be  solved  analytically,  however,  so  the  optimum  expansion  ratio  is  computed 
here  by  finding  the  maximum  impulse  using  a  numerical  search.  The  optimum  expansion  ratio  is 
plotted  as  a  function  of  <j)0  and  r  in  Fig.  7,  and  the  dimensionless  impulse  produced  at  the  opti¬ 
mum  expansion  ratio  is  plotted  in  Fig.  8.  Because  r  is  variable  in  these  figures,  T  is  also,  so  the 
dimensionless  impulse  and  the  average  thrust  coefficient  curves  are  no  longer  directly  propor¬ 
tional  as  they  were  in  earlier  figures,  where  r  was  fixed.  The  optimum  thrust  coefficient  is  plotted 
in  Fig.  9. 

3.3  Variable  Nozzles 

Because  the  expansion  ratio  is  always  optimized  for  variable  nozzles,  the  blow  down  of  vari¬ 
able  nozzles  proceeds  along  the  curve  of  maximum  thrust  in  Fig.  1.  The  dimensionless  impulse 
produced  can  be  computed  as  a  function  of  <j)0  and  r  by  integrating  (11),  and  compared  with  the 

performance  of  fixed  nozzles  operating  at  the  optimum  expansion  ratio.  The  results  are  not  easily 
viewed  in  the  form  of  Fig.  8,  however,  so  they  are  plotted  instead  here  in  Fig.  10  as  the  impulse 
penalty  in  using  optimized  fixed  nozzles  compared  to  the  variable  nozzle  case,  based  on  a  per¬ 
centage  of  the  variable  nozzle  impulse.  The.  penalty  is  seen  to  be  very  modest,  not  exceeding 
about  3%  for  all  cases  calculated.  The  reason  the  penalty  is  not  more  severe  can  be  traced  to  the 
flatness  of  the  curves  in  Fig.  1  near  the  points  of  maximum  thrust.  This  is  illustrated  more  clearly 
in  Fig.  2,  where  the  difference  between  the  variable  and  optimized  fixed  nozzle  blow  down  paths 
is  shown.  As  can  be  seen,  the  maximum  instantaneous  thrust  coefficient  is  not  that  much  higher 
for  the  variable  nozzle  case  than  for  optimized  fixed  nozzles,  leading  to  calculated  average  thrust 
coefficients  which  are  very  close. 

4.0  CONCLUSIONS 

Analytical  solutions  of  the  constant  volume  limit  of  pulsed  propulsion  for  fixed  expansion  ratios 
have  been  explored  and  compared  with  the  case  where  a  variable  nozzle  is  used  to  match  the 
pressure  ratio  at  all  times  during  the  blow  down.  The  solutions  apply  for  all  supersonic  flows 
where  compression  waves  remain  exterior  to  the  nozzle.  It  is  found  that: 

1.  CV  devices  should  optimize  at  the  approximately  the  same  mixture  ratio  as  CP  devices, 
namely  the  mixture  ratio  which  maximizes  the  initial  temperature  and  minimizes  the  mo¬ 
lecular  weight. 


•nU.  **.,  ■  - - 

2.  The  thrust  of  CV  devices  is  maximized  in  the  same  way  as  for  CP  devices,  namely  by 
maximizing  the  initial  pressure  and  the  throat  area. 

3.  The  blow  down  time  depends  on  r  =  pi  p0,  only  weakly  on  y,  and  is  independent  of  the 
pressure  ratio  0O  =  P0  /  P_ . 

4.  In  general  an  optimum  expansion  ratio  exists  which  maximizes  the  impulse  produced  by 
the  blow  down  of  a  CV  device  at  a  fixed  expansion  ratio. 

5.  Using  an  optimized  fixed  expansion  ratio  nozzle  results  in  only  a  modest  impulse  penalty 
compared  to  using  the  more  complicated  variable  expansion  ratio  nozzle.  The  penalty  does 
not  exceed  3%  of  the  impulse  of  the  variable  nozzle  for  all  cases  considered. 
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Fig.  3.  Blow  down  to  r = 0.75 
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Fig.  5.  Blow  down  to  r  =  0.25 


Fig.  6.  Blow  down  times 
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Fig.  7.  Optimum  expansion  ratios 


Fig.  8.  Optimum  impulse 


